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Abstract 

In this paper, starting from some general and plausible assumptions based on ge- 
ometrical optics and on a common feature of the truncated Bessel beams, a heuristic 
derivation is presented of very simple analytical expressions, capable of describing 
the longitudinal (on-axis) evolution of axially-symmetric nondiffracting pulses when 
truncated by finite apertures. We apply our analytical formulation to several situ- 
ations involving subluminal, luminal or superluminal localized pulses, and compare 
the results with those obtained by numerical simulations of the Rayleigh-Sommerfeld 
diffraction integrals. The results are in excellent agreement. The present approach 
can be rather useful, because it yields, in general, closed-form expressions, avoiding 
the need of time-consuming numerical simulations; and also because such expres- 
sions provide a powerful tool for exploring several important properties of the trun- 
cated localized pulses, as their depth of fields, the longitudinal pulse behavior, the 
decaying rates, and so on. 
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1 Introduction 



Ideal nondiffracting pulses (INP) are infinite energy solutions of the ordinary linear 
wave equation, capable of maintaining their spatial shapes indefinitely (sometimes with 
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just small local variations) while propagating [1- 9] . When these ideal solutions are to 
adapted to real situations and applications, they must be spatially truncated by a finite 
aperture (i.e., generated by a finite aperture), getting transformed into finite energy so- 
lutions, with finite field depths, even if these field-depths are very large when compared 
to those of ordinary pulses. 

When we truncate an INP, the resulting wave field cannot be obtained, in general, 
in analytical form. In this case one has to resort to the diffration theory and perform 
numerical simulations of the diffraction integrals, such as that, well known, of Rayleigh- 
Sommerfeld[10]. 

Indeed, one can get very important pieces of information about a truncated 
nondiffracting pulse (TNP) by performimg numerical simulations of its longitudinal 
evolution[ll-14], especially when the pulse is axially symmetric. 

In this paper we shall show that, by using some general and plausible assumptions, 
based on geometrical optics and on a common feature of truncated Bessel beams, a heuris- 
tic derivation of simple analytical expressions is possible, capable of furnishing the lon- 
gitudinal (on-axis) evolution of the TNPs. It is interesting to notice that this approach 
depends only on the spectral structure of the relevant INP. 

We compare the results of our analytical method, when applied to several different 
situations involving subluminal, luminal or superluminal TNPs, with the results obtained 
from the usual numerical simulation of the Rayleigh-Sommerfeld integrals: the results are 
in an excellent agreement. 

This method, due to its extreme simplicity and analytical character, can be a powerful 
tool for exploring several important properties of the TNPs, as their depth of fields, 
the longitudinal pulse behavior, the decaying rates, etc.; for revealing the effects of the 
spectral parameters on the pulses evolution; and also for comparing the "effectiveness" of 
the different kinds of TNPs, as the subluminal, luminal and superluminal ones. Without 
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this method, all those results could be reached (in each particular situation) only by 
performing several time-consuming numerical simulations. 

2 Heuristic approach for describing the on-axis evo- 
lution of TNPs 

Let us begin this Section by making some comments about the truncated Bessel beams 
and about some approximations, which will be used below for developing the method here 
proposed. 

2.1 Some observations and approximations about truncated 
Bessel beams 

An ideal (infinite energy) Bessel beam (IBb) is given by[15] 

*m b (j>,z,t) = J (k p p)e if3z e-^ , (1) 

where p 2 = x 2 + y 2 is the transverse coordinate, k 2 = uj 2 /c 2 — (3 2 is the transverse 
wavenumber, (3 is the longitudinal wavenumber, and u is the angular frequency. 

An important parameter[16-19] of an IBb is its axicon angle 9, where u = c/3/ cos9. 

When a Bessel beam is truncated by a finite aperture of radius R, Eq.Q cannot 
be used to describe the resulting beam in the whole space. However, if the size of the 
aperture is big enough to contain several bright rings of the ideal incident Bessel beam, i.e., 
if R >> 1/kp, we can use geometrical optics to get some characteristics of the truncated 
Bessel beam (TBb) evolution. 

In this case (see Fig.l), we can say that, in the spatial region localized inside the cone 
oriented along the z-axis, with apex at z = Z = Rj tan 9 and base given by the circular 
aperture, the resulting TBb can be approximately described by Eq.(pQ). 
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Figure 1: A typical Bessel beam truncated by a finite aperture. 



However, when using geometrical optics, after the distance Z = Rj tan# the on-axis 
amplitude of the TBb becomes approximately zero (see Fig.l). The distance Z is called 
the depth-field of the TBb. 

Keeping the above observations in mind, one may affirm that, since R >> l/k p , the 
on-axis behavior of a TBb can be approximately described by 



^TBb(p = 0,,M) 

which can be compactly written as 



for z < 



for z > 



tan f 



R 
tan t 



(2) 



TBb 



(p = 0,z,t) 



~ e il3z e~ iu)t 



H(z)-H z- 



R 



tan 9 



(3) 



where H(.) is the Heaviside step function. The equation above is the starting point for 
our heuristic method describing the on-axis TNP's behavior. 

According to the approximation Eq.Q, the on-axis field intensity of a TBb is a 
rectangular function with unitary value, till z = Z. 
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On the other hand, when a numerical simulation of the diffraction integrals is per- 
formed^,^], one can observe that the TBb presents some on-axis field oscillations around 
the unitary value, before suffering an abrupt decay after z = Z . Such oscillations cannot 
be predict by geometrical optics, and arise only due to the abrupt truncation made by 
the aperture. However, it is important to stress here that, despite the fact that those 
oscillations are not predicted by Eq.Q, such an error is not present, in general, in the 
case of our description of truncated localized pulses. One can understand this by noticing 
that, since nondiffracting pulses are constructed through Bessel beams superpositions, 
those oscillations, originating from each TBb, suffer a destructive interference. 

With all this in mind, we are ready to develop our method. 

2.2 The heuristic approach 

It is well known[l-9] that axially-symmetric ideal nondiffracting pulses (INP) can be 
made by zero order Bessel beam superpositions, 

roc r lJ / c / \ ufl \ 

Hp,z,t) = J Q dco J dpS(u t P)J lp\J^-P 2 Je^e-*"* (4) 

provided that the spectral function S(u,(3) entails a linear relationship of the type 

u = V(3 + b (5) 

between ui and (3. In this way, by putting S(u, (3) = S(u))6(u — Vj3 — b), the general form 
of an axially-symmetric INP gets written as 
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where ( = z — Vt, and S(u) is the frequency spectrum. Obviously, the INP will be 
subluminal, luminal or supeluminal, depending on the value of V, it being be <, = or 
> c. The positive quantities u m i n and u) max are the minimum and maximum angular 
frequency allowed for the Bessel beams in the superposition ([B]). and their values have to 
be estimated as follows. 

Once we have chosen the value of V in ©, the values of b, u) min and 0J max are to 
related between them in such a way that 



for all positive angular frequency uj m i n < uj < oj max used in the superposition ©• 

In relation the condition on k p eliminates any unphysical behaviors of Bessel 

functions and evanescent waves in Eq.©. The second condition ((3 > 0) eliminates any 

backwards-travelling Bessel beams in the same superposition, since we are considering 

positive values of the angular frequency only. 
Taking into account conditions ((Zj), one has: 

• For Subluminal (V < c) INPs: b > 0, 0J min = b and 0J max = cbj (c — V) 

• For Luminal (V = c) INPs: b > 0, uj min = b and uj max = oo 

• For Superluminal (V > c) INPs: b > 0, u) min = b and 0J max = oo. Or b < 0, 
Wmin = cbj (c - V) and cu max = oo 

The INPs provided by Eq.© can propagate without distortion indefinitely, with 
peak-velocity V. 

Such INPs possess an infinite energy, and so, for real applications, they must be 
spatially truncated (i.e., generated by finite apertures) [11-14], resulting in finite energy 
solutions, with a finite depth of field. 





/3 > 
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When such truncation is made, the resulting pulse in general cannot be obtained in 
an analytically form, but has to be numerically calculated from the diffraction theory, 
by using, for example, the Rayleigh-Sommerfeld formula[10-14]. That is, once we have a 
known INP solution ^inp, its truncated version ^/tnp, generated by a finite aperture of 
radius R on the plane z = 0, results given by 

r 2lT r R 1 f (z — z') (z-z')) 

* twp( „,m) = I ml w — ( m L_J + [«**]i- r t| (8) 

the quantities enclosed by the square brackets being evaluated at the retarded time 
ct' = ct — D. The distance D = \J {z — z') 2 + p 2 + p' 2 — 2pp' cos (cf> — </>') is the separation 
between the source and observation points. 

Due to its complexity, Eq.(JEJ) has to be solved numerically in most cases. 

Of particular interest is the on-axis behavior of ^tnp- Actually, many important 
information can be extracted from the evolution of ^tnp{p = 0, z, t), like its depth of 
field, the pulse decaying rate, its longitudinal behavior, the effects of the different spectral 
parameters on the pulse evolution; and, even more important, quantity ^tnp(p = 0, z,t) 
can be used to compare the performance of different kinds of TNP, as for example, the 
luminal and the superluminal ones. 

On considering axially-symmetric TNP, and making p = in Eq.(jHJ, we get some 
simplifications, because the integration on <$' can be immediately done: But, even in this 
case, the integration on p' rarely can be carried on analytically, due to the complexity of 
the integrand, and numerical simulations are once more required. 

To overcome this problem, let us propose the following heuristic approach: 

First, we consider the Bessel beam superposition (JHJ), which provides us with the 
INPs. 

Second, we make the assumption that each Bessel beam, with frequency u and axicon 
angle 9(u), entering in superposition (0), obeys the following condition: 
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kp usm9(u) ' 

R being the radius of the finite aperture that will be used for truncating the INP. 

The above assumption is very plausible, since efficient TNPs are generate by large 
apertures. 

Once condition Q is fulfilled by all Bessel beams in superposition (JHJ), we can use 
again the geometrical optics and assume that, after the truncation, the on-axis behavior 
of each one of those Bessel beams can be approximated by Eq.fjSJ). 

Third, taking into account Eqs. ((HJ) and we may conjecture that the on-axis evo- 
lution of the truncate nondiffracting pulse is approximately given by: 



^tnp(p = 0,z,t) ~ e 



-ibz/V 



duj S(uj) e 



iui(/V 



H(z) -H\z 



R 



, (10) 



tan 9{uj) 

where H(.) is the Heaviside step function and, let us recall, 9{uj) is the axicon angle of 
the Bessel beam with angular frequency u. 

We should notice that in the integrand of (JTUJ), the step function H(z — Rj tan^o;)) 
depends on uj, through of #(u;). 

We can rewrite Eq.()l(JJ) in the form 



q TNP (p = 0,z>0,t)~ e 



-ibz/V 



S{u) e iuJ<:/v du 



/Umax . , 

S(u) e l ^ /v H 
— i m i n 



R 



tan6 l (a;) 



duo 



where the first term in the r.h.s. of (JTTjl is nothing but the INP ^inp{p = 0,z,£), while 
the second term is the perturbation due to the truncation. 
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Now, remembering that for a Bessel beam of axicon angle 6 we have u = cf3 / cos 6, and 
that the spectra of INPs impose the constraint uo = V (3 + b between angular frequencies 
and longitudinal wavenumbers, it becomes easy to show that 



and thus 



R 



tan#(d;) 



\ 



R 



1 - ( — — — 



= - 

2 XV 



bc_" 
Vu 



(12) 



Hz 



R 



tan 9(uj) 



H z 



1 for z > 



for z < 



i- 



\/ 1- (F—fe) 2 



(13) 



With all what precedes, we can eventually write 



^ttvp(p = 0,z,t) ~ e" 



-ibz/V 



5(cj)e iwC/y ^ 



I 



\ 



c_ bc_\ \ 

V VujJ 



V VujJ ) 



du 



(14) 



In the next subsection, we will analyze the fundamental Eq. (jl4}) for the three possible 
types of TNPs: subluminal, luminal, and superluminal. 
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2.2.1 Subluminal TNP 



For the subluminal pulses (V < c), we have b > 0, 0J min = b and 0J max = cbj (c — V). 
In this way, taking into account these facts, and that z > and w min < u < u max , 
we can show, after several manipulations, that Eq. (jl3j) can be written as 



H 



z — 



I c_ be 

\V Vuj 



1 for uj < 



be 



zV 



c — 



for uj > 



Vz 2 + R 2 
be 



zV 



Vz 2 + R 2 



(15) 



Now, by noting that uj min = b < be/ (c — zV j \/ z 2 + R?) < u) max = be/ (c — V), one can 
write Eq.fJHJ), for the subluminal 



y TNP (p = 0,z>0,t) ~e 



-ibz/V 



be 
c-V 



S(u) e^ /y du - f c VOT S{uj)e iujC/v du 

J b 



-ibz/V 



S{u) e l ^ /v du 



\/z 2 +R 2 



(16) 



which represents our method in the case of subluminal TNPs. It is a very 
simple equation, capable of providing closed-form, analytical results for several different 
frequency spectra S(u), as we shall see in Section 3. 

2.2.2 Luminal TNP 



For luminal TNPs (V = c), we have b > 0, 0J min = b and 0J max = oo. 
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With this, and taking into account that z > and LU min < to < uj max , we can, after 
several manipulations, show that Eq. (jl3j) may be written as 



H 



z — 



i I c_ be 



1 for lu < 



1 - 



for uj > 



Vz 2 + R 2 
b 



(17) 



1 - 



Vz 2 + R 2 



and, by noting that uj min = b < b/(l — zj\J z 2 + R 2 ) < u max = oo, we can write Eq. (fT4l 
for the luminal 



q TNP (p = 0,z>0,t) ~e 



-ibz/c 



S(uo)^ /c duo - 1 S{uj)^ /c duj 



-ibz/c 



S(u)e iuC / c du 



\/z 2 + R 2 



Equation (|18j) . which represents our method in the case of luminal TNPs, 
is very simple too, and can provide closed-form, analytical results for many different 
frequency spectra S (to), as we shall see in Section 3. 



2.2.3 Superluminal TNP 

For superluminal TNPs (V > c), the value of b, in the the spectral constraint (0), can 
assume negative or positive values, i.e., — oo < b < oo. Let us analyze the superluminal 
case of Eqs. (fT3|) and (fT^) for both situations, b < and b > 0. 

Superluminal case for b < 
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In this case, we have ui min = cbj (c — V) and u) max = oo. 

Taking into account that z > and 0J min < u> < u> max , and, again, after several 
manipulations, we can show that for this case Eq. ljlHj) can be written as 



H 



r? ( c be 



1 " (v ~ vt) J 



1 for u) < 



be 



zV 



1 for uj > 



Vz 2 + R 2 
be 



zV 



and z < 



and z > 



c — 



Vz 2 + R 2 



R 

V 2 /c 2 - 1 
R 

V 2 /c 2 - 1 



otherwise 



(19) 



By noting that, when z < R/JV 2 /c 2 — 1, we have bejic — zV/\/z 2 + R 2 ) < and 



that, when z > R/\jV 2 /c 2 — 1, we have bc/(c — zV/ \J z 2 + R 2 ) > uj min = bejic — V), one 
can write Eq. (fH|) . for the case V > c and b < 0, as 



^tnp{p = 0,z>0,t) 



roo 

e bc S{u)e iul V v du 



for z < 



v /y2 /c 2_ 1 



/ C e iuC/V doJ for z y 



bc 
c-V 



v /y2/ c 2_ 1 



(20) 



Equation (fTSj) . represents our method in the case of superluminal TNPs with 
b < 0. Again, the integrals are very simple, and can provide closed-form, analytical results 
for many spectra S(u). 
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Before going on to the next case, one can immediately see from Eq. (J2(J|) that, 
independently of S (a>), the superluminal TNPs with b < will reach the distance 
z = Rj JV 2 /c 2 — 1 without deforming. 



Superluminal case for b > 

In this case, we have u min = b and u max = oo. 

Remebering that z > and u m i n < to < uj max , and, as before, after several ma- 
nipulations, we can show that in this case Eq.([13|) can be written in the same form of 
Eq.(HnD. 

Now, taking into account that, when z < Rj ^Jv 2 /c 2 — 1, it is 0J min = b < 
bc/(c — zV/^/z 2 + R?) < ijJ m ax = oo and that, when z > Rj \jV 2 /c 2 — 1, one has 
be/ (c — zV I V z 2 + R?) < 0, we can write our fundamental Eq. (fPf|) . for the case V > c and 
b > 0, in the form 



q TNP (p = 0,z>0,t) 



-ibz/V 



/oo 
S{u)e l ^ /v duo for z< 



R 



V / V 2/ C 2_ 1 



\/z 2 +R 2 



for z > 



R 



v /y2/ c 2_ 1 



:2D 



Equation (|18jl. represents our method in the case of superluminal TNPs with 
b > 0. Again, many closed-form results, for many different spectra S(uj), can be obtained 
from the above equation. 

We can also notice from Eq.()21|) that, for V > c and b > 0, the superluminal TNPs 
will get very low intensities after the distance z = Rj ^V 2 /c 2 — 1. 

It is important to notice that in our method, given by Eqs.(fH)J). (fTHj) . (J2UH2TJ1 . the on- 
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axis evolution of a TNP depends only on the frequency spectrum S(uj) of its corresponding 
INP ^inp, at variance with the Rayleigh-Sommerfeld formula (JHJ), which depends on the 
mathematical expression of ^inp- 

Now, we shall go on to the next Section, where our method will be applied to some 
important kinds of localized waves, and our results will be compared with those obtained 
by making numerical simulations by the Rayleigh-Sommerfeld formula. 

3 Application to some important cases. Closed-form, 
analytical results, and their comparison with nu- 
merical simulations of the Rayleigh-Sommerfeld 
formula. 

The method we have developed in the previous Section is described by Eqs.(|16|). 
(fT£J) . (ffiH 121]) . for the cases of truncated subluminal, luminal and superluminal pulses, 
respectively. 

We are going to apply, now, this method to some important cases involving TNPs; and 
it will be shown that the results agree very well with those obtained through numerical 
simulations of the Rayleigh-Sommerfeld formula (|5jl. 

3.1 The method applied to a truncated subluminal pulse 

A well known ideal subluminal (V < c) nondiffracting pulse [20] is 




where 7 = l/J 1 — V 2 /c 2 , rj = z — c 2 t/V and, as always, ( = z — Vt. 

Now, we want to describe the on-axis behavior of the truncated version (^tnp) of the 
the ideal solution (J22J)- 
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The subluminal INP (I22j) is generated by the superposition with a constant 
spectrum 5" (a;) = c/(25V7 2 ). However, we have to notice that this solution possesses 
backwards-travelling components: Actually, it has uj m i n = bc/(c + V), instead of uj min = b 
(which assures forward components only in ©, as we have seen in Section 2). It is the 
price paid to get such a closed-form and exact INP solution. 

Anyway, we may, and we must, minimize the contribution of those "backwards" com- 
ponents by choosing the subluminal velocity V in such a way that (c + V)/(c — V) » 1. 
Once this condition is satisfied, we can observe that the INP ()22|) is, then, similar to that 
which would be obtained with the same S(u>), but with uj m in = b. 

It should be noticed that the comments and observations above have nothing to do 
with our method, which is constructed since the beginning in order to comprehend causal 
(forward) solutions only, and can be used for any values of the velocity V. Those remarks 
were made just in order that a causal behavior of the INP (j22j) be guaranteed. 

Now, on using S(u) = cj (26V^7 2 ) in our Eq.()16|). which describes the on-axis evolution 
of the subluminal TNPs, one gets 



TNP 



(p = 0,z> 0,t) 



2bVj 2 



-ibz/V 



^IV duJ 



\/z 2 + R 2 



cV 
2bV 1 2 iC 



-ibz/V 



exp 



ibc 
V{c-V) 



C - exp 



ibcVz 2 + R 2 



V{c\/z 2 + R 2 - zV) 
(23) 

which is a very simple closed-form, analytical expression. 

First, we use Eq. fl2"3*j) to get the pulse-peak intensity behavior. To do this, we just 
put z = Vt -»• C = in (gBD - 

Let us consider three different cases: (1) V = 0.995c and b = 1.5 x 10 15 Hz; (2) 
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V = 0.998c and b = 6 x 10 14 Hz, and (3) V = 0.9992c and b = 2.4 x 10 14 Hz. In all 
cases, we consider the radius of the finite aperture to be R = 4 mm. 

At the same time, we compare our result with that obtained by the numerical simu- 
lation of the Rayleigh-Sommerfeld formula (JEJ), with ^inp given by ([22)1 . 

Figure 2 shows the plots, where those represented by continuous lines have been 
obtained from our equation (|23jh while those represented by dotted lines come from the 
numerical simulation of (JEJ). 

We can verify the excellent agreement existing among those results. 




0.1 0.15 

z [m] 



0.2 



0.25 



Figure 2: The peak intensity evolution of this subluminal TNP for the three cases: (1) 
V = 0.995c and b = 1.5 x 10 15 Hz. (2) V = 0.998c and b = 6 x 10 14 Hz. (3) V = 0.9992c 
and b = 2.4 x 10 14 Hz. In all cases R = 4 mm. The continuous lines are obtained from our 
closed-form analytical expression (J23J) and those represented by dotted lines come from 
the numerical simulation of the Rayleigh-Sommerfeld formula (JHJ). 
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Now, we are interested in the on-axis longitudinal pulse evolution, in the three cases 
above considered, in the time instants given by t' = 0.11 ns, t" = 0.22 ns and t'" = 0.33 ns. 

Figures (3a,3b,3c) show the results corresponding to the cases (1), (2) and (3), respec- 
tively. As before, the continuous lines represent the results obtained from our Eq. (j23|) . and 
the dotted ones those coming from the numerical simulation of (|8jl . Again, we consider 
R = 4 mm. 




Figure 3: The on-axis evolution of this subluminal TNP, in the times t' = 0.11ns, t" = 
0.22ns and if" = 0.33ns, for each case cited in Fig. 2. Figure 3a, 3b and 3c represent the 
cases (1),(2) and (3) respectively. The continuous lines are the results obtained from our 
closed-form analytical expression ()23|) and those represented by dotted lines come from 
the numerical simulation of the Rayleigh-Sommerfeld formula (jSJ). 

We can observe, once more, a very good agreement among the results, confirming 
that our method works very well. 

3.2 The method applied to the truncated luminal Focus Wave 
Mode pulse 

A very famous ideal luminal (V = c) nondiffr acting pulse[l] is the Focus Wave Mode 
pulse (FWM), given by 
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VlNp(p,Z,t) 



ac — i( 



ac 



exp 



( 



-ib 



2c 



7] exp 



( 



2c[ac — i() 



-bp 2 



(24) 



where r] = z + ct, ( = z — ct, and a > is a constant. 

Like all the INPs, the FWM possesses infinite energy and must be truncated (i.e., 
generated by a finite aperture) for real applications. We shall use our method to get closed- 
form, analytical expressions for the on-axis evolution of the truncated version ^tnp of 
Eq.flMI). 

The exact ideal solution is obtained from the superposition (JHJ) with V = c 
and S(lj) = aexp(a6/2)exp(— aw). Because it has u) min = 6/2, instead of u min = b, its 
spectrum possesses backwards components [1-6] in the range b/2 < u < b. To overcome 
this problem, we must minimize the contribution of the nonphysical part of the spectrum; 
and this can be done if ab << 1. Once such a condition is obeyed, the INP ()24|) can be 
considered similar to the one that we would obtain with the same frequency spectrum of 
the FWM, but with uj min = b. 

Again, let us notice that such remarks are made just to validate the causality of the 
INP ()24j1 . and have nothing to do with our own method. 

Now, by using 5"(a;) = aexp(a6/2)exp(— aw) in our Eq.ffTHj). we get the on-axis evolu- 
tion of the truncated FWM: 



which is a very simple closed-form, analytical expression. 

First, let us put ( = in our Eq. (j25|) in order to get the pulse-peak intensity behavior. 




(25) 
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We consider three different cases: (1) a = 1.6 x 10 16 and b = 5 x 10 11 Hz; (2) 
a = 1.25 x l(T 16 s and b = 3 x 10 n Hz; and (3) a = 1 x 1(T 16 and b = 2 x 10 n Hz. In all 
cases we adopt the aperture radius R = 2 mm. 




z [m] 

Figure 4: The peak-intensity evolution of the truncated luminal FWM pulse for the three 
cases: (1) a = 1.6 x 1(T 16 and b = 5 x 10 n Hz; (2) a = 1.25 x l(T 16 s and b = 3 x 10 n Hz; 
(3) a = 1 x 10~ 16 and b = 2 x 10 11 Hz. In all cases R = 2 mm. The continuous lines are 
obtained from our closed-form analytical expression (J23j) and those represented by dotted 
lines come from the numerical simulation of the Rayleigh-Sommerfeld formula (JSJ)- 

Figure 4 shows the results. The continuous lines represent the results obtained from 
our Eq. (j25|) . while the dotted ones are the results of the numerical simulation of the 
Rayleigh-Sommerfeld formula (|S|). 

The results agree so well, that the corresponding continuous and dotted curves do 
superpose to each other. 

Now, we are going to use our Eq. (|25|) to show the on-axix evolution of this TNP, in 
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the three cases above considered, for the time instants t' = 0.22 ns, t" = 0.44 ns and 
t'" = 0.66 ns. 

Figures (5a,5b,5c) show the results corresponding to the cases (1), (2) and (3), re- 
spectively. The continuous lines are the results obtained from our Eq.(j2SJ), and the dotted 
lines those coming from the numerical simulation of (JHJ). Again, we consider R = 2 mm. 




z-ct[m] slO" 7 z-ct[m] xlO" 7 z-ct[m] x 10" 7 

(a) (b) (c) 



Figure 5: The on-axis evolution of the truncated luminal FWM pulse, at the times i! = 
0.22 ns, t" = 0.44 ns and t'" = 0.66 ns, for each one of the cases cited in Fig.4. Figures 
5a, 5b and 5c represent the cases (1),(2) and (3), respectively. The continuous lines are the 
results obtained from our closed- form analytical expression (J23|) . while those represented 
by dotted lines come from the numerical simulation of the Rayleigh-Sommerfeld formula 

The results are in excellent agreement, showing the very good efficiency of the method. 

3.3 The method applied to the truncated Superluminal Focus 
Wave Mode pulse 

An interesting, approximated, superluminal (V > c) ideal nondiffr acting solution to 
the wave equation is the so-called[6] Superluminal Focus Wave Mode pulse (SFWM): 



9 mP (p,z,t) = aVexp(^) X exp ( ^±^_ (( a V - <) - X' 1 )) , (26) 
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where rj = z + Vt, ( = z — Vt, a > is a constant, and where 

X = [iaV-i0 2 +{^-l)p 2 ) ' (27) 

Expression (j21j) is a very good aproximate solution of the wave equation[6] if ab « 1, 
which is also the condition for minimizing the contribution of the backwards components 
of pijl. Actually, this superluminal INP can be obtained from superposition ©, with 
b > 0, when using 5(o;) = aexp(a6/2)exp(— au), with constant a > 0, but with u min = b/2 
instead of u) min = b. 

To get the closed-form, analytical mathematical expression that describes the on-axis 
evolution of the truncated version of PS)), let us put S(u) = aexp(ab/2)exp(—auj) in our 
Eq.(j2H): 



^tnp(p = 0,z>0,t) 



a y e ab/2 e ~ibz/V 

aV-iC 



exp 



-bcVz 2 + R 2 (aV-iCY 
VicV^TK 1 - zV) 



for z < 



R 



^V 2 /c 2 -l 



for z > —, — = , 



(2f 



Now, let us set C = in our Eq.(|28|). for analyzing the peak-intensity behavior of the 
truncated SFWM. 

We consider three different cases: (1) V = 1.0002 c, b = 3 x 10 12 Hz and a = 2.5 x 
10~ 17 s; (2) V = 1.0001 c, b = 1 x 10 12 Hz and a = 5 x 10~ 17 s; and (3) V = 1.00008 c, 
b = 2 x 10 12 Hz and a = 1.1 x 10~ 17 s. In all these cases we choose R = 3 mm as being the 
radius of the aperture. 

The plots are shown in Fig. 6, where the continuous lines represent our results using 
Eq.(|28|). and the dotted ones represent those obtained from the numerical simulation of 

21 



the Rayleigh-Sommerfeld formula (jSJ). We can observe an excellent agreement among the 
results. 

l^(p=o,z = vt)| 2 




z [m] 

Figure 6: The peak intensity evolution of the truncated superluminal Focus Wave Mode 
pulse for the three cases: (1) V = 1.0002 c, b = 3 x 10 12 Hz and a = 2.5 x 10~ 17 s; (2) 
V = 1.0001 c, b = 1 x 10 12 Hz and a = 5 x 10~ 17 s; and (3) V = 1.00008 c, b = 2 x 10 12 Hz 
and a = 1.1 x 10~ 17 s. In all cases R = 3 mm. The continuous lines are obtained from 
our closed- form, analytical expression (J28j) . while those represented by dotted lines come 
from the numerical simulation of the Rayleigh-Sommerfeld formula (JHJ) 

Now, we are going to use our method to show the on-axis evolution of the pulse 
intensity at three different times t' = 0.14 ns, t" = 0.29 ns and t'" = 0.43 ns, for each one 
of the cases cited above. 

Figures (7a,7b,7c) show the plots. Again, the curves given by continuous lines come 
from our Eq.([28jh while the dotted ones from the numerical simulation of (|8*jl. 

Once more, there is an excellent agreement among the results, confirming the validity 
and efficiency of our method. 
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Figure 7: The on-axis evolution of the truncated SFWM pulse, at the times t' = 0.14 ns, 
t" = 0.29 ns and t'" = 0.43 ns, for each one of the cases cited in Fig.6. Figures 6a, 
6b and 6c represent the cases (1),(2) and (3), respectively. The continuous lines are the 
results obtained from our closed- form, analytical expression (|25jl. and those represented 
by dotted lines come from the numerical simulation of the Rayleigh-Sommerfeld formula 



4 Extending the present method to the almost non- 
diffracting (finite energy) pulses, truncated by fi- 
nite apertures 

In the previous Sections we have developed a (heuristic) method capable of providing 
closed-form analytical expressions, describing the on-axis evolution of the INP truncated 
by finite apertures. 

It is well known [1,6] that, besides the INPs, there are the almost nondiffr acting pulses 
(ANP), which also need infinite apertures to be generated, but possess a finite energy 
content. 

Once a function S(u) is chosen, and an INP with a velocity V and b = bo is obtained 
from Eq.©, we can get an ANP by integrating Eq.(|B*)l. over the parameter b, with a 
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suitable choice of the weight function S'(b) which has to be concentrated around b = b . 
More explicitly: 



(29) 

where S"(u, b) = S(u)S'(b) is a spectral function with S"(6) well localized around b = b . 
Obviously, one can recover the INPs just by adopting the choice S"(u, b) = S(u)5(b — 

bo). 

An ANP can be viewed as a Bessel beam superposition (0J) with a spectral function 
S((j,(3) well concentrated around a straight line uj = V (3 + bo. 

The ANPs are interesting solutions, due to their finite energy contents, and can 
maintain their spatial shape for long distances[l,6]. 

However, even possessing finite energy, the ANPs need — as we were saying — infinite 
apertures to be generated: something that cannot be obtained in the real world. Due to 
this fact, it is rather important to know the behavior of these pulses when truncated by 
finite apertures; that is, to know the ^tnps versions of the ^anps- This can be got by 
making numerical simulations, again, of the Rayleigh-Sommerfeld integral formula (jHj). 
on replacing ^inp with anp- 

On the other hand, the extension of our method to the cases of ANPs 
truncated by finite apertures can be performed in a very simple way, just by 
multiplying our fundamental equations (|16|) . (|18|) and (|2()l I21|) by the S(b') un- 
der consideration, and performing the relevant integration over the parameter 
b. 

We are going to show this by the following example, in which we shall obtain closed- 
form analytical expressions for the truncated version ^tnp with a well known finite energy 
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ANP[1,20]: Namely, the Modified Power Spectrum pulse (MPS). 



Example: The truncated version of the MPS pulse. 

A well-known luminal ANP is the MPS pulse, which can be obtained by integrating, 
over the parameter b, the FWM pulse (J23|) with the weight function 

S'(b) = H(b - b ) qexp(-q(b - b )) , (30) 

quantities q and b being positive constants. More explicitly, the MPS pulse can be written 
as: 

= £ ^ (^) ^ Uac -zC) ) eM ~ q{b ~ bo)) dh 

(31) 

= 2ac 2 q f -ibp \ ( b p 2 \ 

~ (2cq + irj)(ac - i() + p 2 6XP \ 2c V J 6XP \ 2c(ac - i() ) ' 

with rj = z + ct, C = z — ct. 

This ANP has a finite energy content; however, it needs an infinite aperture to be 
generated. 

We shall use the extended version of our method to get a closed-form analytical 
expression for the on-axis evolution of the truncated version, ^tnp, of the MPS pulse. 
As we have seen, to get this we just need multiplying the truncated version of the FWM 
pulse, given by Eq.()25|). by the corresponding weight function S'(b) given by ()3U|). and 
performing the integration over the parameter b. In this way, the on-axis evolution of the 
truncated MPS pulse is given by: 
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»™,(, = 0, * > 0, t) = r -^e°V V <^exp(- # - W)«p( ~-^™^f 

Jb ac-iQ \ c(\/z 2 + R 2 - z) 



aqce 



ab /2 e -ib {) z/c 



, ( a iz V z 2 + R 2 (ac — i() 

(ac _ lC) ^__ + _ + __ r _ 



exp 



W-2 2 + R 2 {ac-iCY 



iVz 2 + R 2 



(32) 



which is a closed-form analytical expression. 

As before, let us put ( = in our Eq.((H2j) in order to get the pulse peak intensity 
behavior. 

riV(p=o,z=ct)i 2 




Figure 8: The peak-intensity evolution of the truncated luminal MPS pulse for the three 



-16 



, fro 



5 x 10 Hz and q 

-i-iic 



2 x 10- n s; (2) a 

-16 



1.25 x 10 



-16, 



cases: (1) a — 1.6 x 10 

b = 3 x 10 n Hz and q = 10 x 10- n s; and (3) a = 1 x 10~ 16 , b = 2 x 10 n Hz and 
q = 20 x 10~ n s. In all cases R = 2 mm. The continuous lines are obtained from our 
closed- form analytical expression (|32jl . while those represented by dotted lines come from 
the numerical simulation of the Rayleigh-Sommerfeld formula (JEJ). 
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Let us consider three different cases: (1) a = 1.6 x 10 16 s, bo = 5 x 10 11 Hz and 
q = 2 x Kr n s; (2) a = 1.25 x l(T 16 s, & = 3x 10 11 Hz and q = 10 x 10- n s; and (3) 
a = 1 x 10~ 16 s, bo = 2 x 10 11 Hz and q = 20 x 10 -11 s. In all cases, we adopt the aperture 
radius R = 2 mm. 

Figure 8 shows the results. The continuous lines represent those obtained from our 
Eq.(|32|). while the dotted ones are the results of the numerical simulation of the Rayleigh- 
Sommerfeld integral formula (JHJ). 

One can verify the excellent agreement among the results. 

Now, we are going to use our Eq. (|32J1 to investigate the on-axis evolution of this TNP 
in the three cases above considered, for the instants t' = 0.22 ns, t" = 0.44 ns and 
t'" = 0.66 ns. 

Figures (9a,9b,9c) show the results corresponding to the cases (1), (2) and (3), re- 
spectively. The continuous lines come from our Eq. ()32)1 . while the dotted lines are those 
coming from the numerical simulation of (|8*J). Again, we consider R = 2 mm. 
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Figure 9: The on-axis evolution of the truncated luminal MPS pulse, at the times 
t' = 0.22 ns, t" = 0.44 ns and t'" = 0.66 ns, for each one of the three cases consid- 
ered in Fig. 8. Figures 9a, 9b and 9c represent the cases (1),(2) and (3), respectively. 
The continuous lines are the results obtained from our closed-form analytical expression 
(I32J) . while those represented by dotted lines come from the numerical simulation of the 
Rayleigh-Sommerfeld formula (|5|). 
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Once more, there is an excellent agreement among the results, confirming the validity 
and efficiency of our method. 

Before finishing this Section, it is important to notice that the closed-form analytical 
expressions of the truncated ANP, obtained with our method, can be advantageously used 
for comparison purposes with the corresponding non-truncated ANP: thus illustrating the 
effects due to the truncation, and, for example, telling us till what distance we can use 
the 3D solution of the ANP as a good approximation to the corresponding 3D TNP. 

5 Conclusions 

In this paper we have developed a very simple method for describing the space- 
time on-axis evolution of truncated nondiffracting pulses, be they subluminal, luminal or 
superluminal. 

It is important to notice that in our method, given by EqB.(|IBj). (fTHj). and (}27fl |2*T|) . 
the on-axis evolution of a TNP depends only on the frequency spectrum S(uj) of its 
corresponding INP ^jnp, contrary to the Rayleigh-Sommerfeld formula (jHJ) which depends 
on the explicit mathematical expression of ^inp- We also have extended our method to 
describe the truncated versions of the ANPs. 

Due to such a simplicity, we can obtain closed-form analytical expressions, which 
describing the on-axis evolution of innumerable TNPs. In this paper we have done that 
for the truncated versions of a few, very well-known, localized waves: subluminal, luminal, 
or superluminal. 

We have compared our results with those obtained through the numerical simulation 
of the Rayleigh-Sommerfeld integrals, and we have observed an excellent agreement among 
them, confirming the efficiency of our method. 

The present approach can be very useful, because it furnishes, in general, closed-form 
analytical expressions, avoiding the need of time-consuming numerical simulations; and 
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also because such expressions provide a powerful tool for exploring several important 
properties of the truncated localized pulses: as their depth of field, the longitudinal pulse 
behavior, the decaying rates, etc. 
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